Meson loops and the go'D-n coupling 
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The D* Div form factor is evaluated at low and moderate Q 2 in a hadronic loop model, for off-shell 
D mesons. The results contain arbitrary constants, which are fixed so as to match previous QCD 

^> ■ sum rule calculations valid at higher Q 2 . The form factor obtained in this way was used to extract 

the coupling constant, which is in very good agreement with the experimental value. 
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I. INTRODUCTION 



The measurement of the D*Dir coupling, made by the CLEO collaboration [1], yelding the result gwD-n — 17-9 ± 
0.3 ± 1.9, created embarrassement in the community of QCD sum rules (QCDSR) which, by means of various schemes, 
had predicted much smaller values. Several different approaches of QCDSR were employed: two point function 
combined with soft pion techniques [2,3], light cone sum rules(LCSR) [4,5], light cone sum rules including perturbative 
£C) ' corrections [6] , sum rules in a external field [7] , double momentum sum rules [8] and double Borel sum rules [9] . The 
LCSR prediction made in ref. [4] became even smaller after the radiative corrections have been included [6]. The 
upper limit of these predictions was gD'D-n = 13.5 [6] and it turned out to be 30% lower than the central value of the 
CLEO measurement. 

Although the QCD sum rule approach certainly suffers from large uncertainties, in several other cases good agree- 
On] ments with experiment were obtained. Therefore, we cannot simply be skeptical about the whole sum rule approach. 
(JLj ■ The qd'D-k coupling constant does not seem, a priori, to be particularly exotic and other theoretical approaches did 
not produce such large discrepancies for this quantity. A careful discussion meant to reduce the uncertainties presented 
, by quark models, performed in the framework of Dirac equation [10], and prior to the experimental measurement, 
has led to the result gwD-x — 18. It should be stressed that this result has been obtained in the heavy quark limit. 
The recent (quenched) lattice QCD calculation has produced gD*Dn = 18.8 ± 2.3+^0 [H]- R 1S therefore important 
to understand the specific difficulty which the standard sum rule approach seems to encounter in this case. 

After the appearance of experimental data, three works [12-14] tried to reconcile the LCSR estimates with the 
measured figure. In [12], it was noted that the inclusion of an explicit radial excitation contribution to the hadronic 
side of the LCSR (often referred to as the left hand side of the sum rule) could significantly improve the value of 
9d*Dtx and, at the same time, the stability of the sum rule with respect to the Borel parameter M . 

In the standard QCDSR approach of [9], a modification in the continuum contribution (such as the explicit inclusion 
of a radial exictation) does not seem to be neither necessary nor promising, because there the Borel suppression is 
much more effective. In ref. [9] the go* Dn (Q 2 ) form factor was estimated as a function of the off-shell pion momentum 
Q 2 . Since the sum rule obtained from the three-point function adopted is not valid at Q 2 = 0, in order to determine 
the D*Dir coupling, it was necessary to extrapolate the Q 2 behaviour of the form factor. Of course there are large 
uncertainties in this procedure, and, the value of go* d-k = 5.7 ± 0.4 obtained was much smaller than the experimental 
result. In a subsequent calculation of the DDp vertex [15], the DDp form factor was calculated for both D and p off- 
shell mesons and the QCDSR results were parametrized by analytical forms such that the respective extrapolations to 
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the D and p poles provided consistent values for the gDDp coupling constant. This method of double parametrization 
plus matching at the on-shell point was then employed in [16] to recalculate the D*Dir coupling and led to the value 



9D'Dn = 14.0 ± 1.5 . 



(1) 



While this number is much closer to the experimental value, there is still a discrepancy. Moreover, the procedure 
of fitting the QCDSR points in the deep euclidean region and extrapolating them to the time-like region still contains 
uncertainties, such as, for example the analytical form chosen for the parametrizations, i.e., monopole, exponential or 
gaussian. 

In the present work we return to this question and employ hadronic loops, calculated by means of effective field 
theories (EFT), in order to produce a better parametrization for D*Dw results derived by means of QCDSR. Purely 
hadronic calculations are independent from QCDSR and involve the choice of an effective Lagrangian, including the 
possible requirements of chiral symmetry and/or SZ7(4). However, beyond tree level, one has to deal with the problems 
and uncertainties associated with renormalization. As we discuss in the sequence, a suitable combination of EFT and 
QCDSR results allows the elimination of undcsircd indcterminacies of both approaches, improving significantly their 
predictive powers. Effective interactions are discussed in section II, whereas results and conclusions are presented in 
sects. IV and V. As an straightforward exercise we also make a prediction for gB*B-K- 



The full D*Dir vertex function is shown in Fig. 1. Leading contributions to this vertex come from both the tree 
interaction and the three classes of diagrams depicted in Figs. 2, 3 and 4. Meson loops are a necessary consequence 
of quantum field theory and do contribute to several hadronic observables. In practice, due to problems associated 
with infinities, renormalization becomes unavoidable in the evaluation of loop corrections to observables. Nowadays, 
this kind of procedure is rather well established at the hadronic level, in processes such as pion-pion and pion-nucleon 
scatterings. In the case of D mesons, on the other hand, the theory is much less developed and hence we resort to an 
alternative. The basic idea is to isolate the unknown loop parameters into some basic constants, in such a way that 
they can be determined by matching the results of loop and QCDSR results. 

Before proceding, some remarks are in order. We first note that some diagrams, such as, for instance, that in 
Fig. 2a , contain internal vertices involving the D*Dir coupling. This suggests that the calculation is "cyclic", since 
one needs to use the D*Dtt form factor in order to calculate the D*Dtt form factor. Actually, there are differences 
between the internal particles and the external ones. The former are always virtual, whereas the latter may be either 
real or put on mass shell in the extraction of the coupling constant. In the framework of perturbation theory, at 
leading order, internal particles arc treated as elementary, without structure. They are assumed to be point-like and 
the evaluation of leading terms does not require the use of internal form factors. Consistently, one must use bare 
coupling constants for these interactions. 

There are heavy mesons circulating in the loops shown in Figs. 2-4 and one might be tempted to argue that other 
states should also be included. We do have, for example, fermion-antifcrmion components such as NN or A C A C in the 
loops. An incoming positive pion can split into a p plus a n, and so on. However, in a different context [17], it has 
been shown that this kind of splitting is suppressed with respect to the pion — > meson-meson splitting, by one order 
of magnitude. The neglect of this kind of contribution seems therefore justified. The same holds for the possibility of 
strangeness circulating in the loop, associated with virtual states such as D Sl D* sl K and K*. Using only 7r's, p's, D's 
and D*'s we cover the low and high Q 2 regions of the form factor. Thus it is enough to our purposes to work with a 
simple effective theory, involving only 7r's, p's, D's and D*'s that, as has been discussed elsewhere, [18], has proven 
to be phenomcnologically successful. The same happens in this work. 

The diagrams considered in this calculation have been divided into three classes. The first one, represented in Fig. 
2, involves only triangular loops. The processes in Fig. 3, on the other hand, contain bubble-type loops and four-leg 
vertices associated with the gauge structure of p interactions. Finally, diagrams given in Fig. 4 involve bubble loops 
as well as single-particle propagators. This last feature might suggest that these interactions should be considered as 
mass corrections. However, the nature of the effective interactions described below is such that, in some terms, the 
poles of the single-particle states are cancelled. These terms do correspond to proper three-point functions and are 
kept in the evaluation of the form factor. Pole cancellations are indicated by crosses in Fig. 4. 

We adopt an effective Lagrangian constrained by SU (2) flavor and chiral symmetries, as well as gauge invariance. 
The ttDD* interaction is given by [18-21] 



II. EFFECTIVE DYNAMICS 



£->rDD» — i g-rrDD 



.. [Dr a D; - D*T a D] d^ a . 



(2) 
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where r a are the Pauli matrices, <j> a denotes the pion isospin triplet, while D = (D°,D + ) and D* = (D*° , D* + ) 
represent the pseudoscalar and vector charm meson doublets, respectively. The hat on top of the coupling constant 
indicates its bare nature. 

The p couplings are assumed to be universal and are implemented by covariant derivatives of the form 



= d» - ig p T-ff 



(3) 



where g p is the universal coupling constant and T is the isospin matrix suited to the field it is acting upon. In 
this work we need 



V^D = D 



D 



Wfa = [d»+g p e abc p^} cf> c . 
Using this prescription in cq. (2), we obtain 

£pnDD* = i 9p ()nDD* [DT b D*^ - D' p T b D] E bca p£<j> a 

For the other couplings, we depart from the free Lagrangians and have 

Attt = \d p (p a d^(j)a ~> £ p7 ™ = 9p € abc 4> a dp4>b Pc i 



C DD = d p D d^D - C pDD = - g p [D t c d p D - 8 P D r c D] p£ , 



£ D * D * 



-pD*D* 



-- [d p E>l - ft,D*] [d»D* v - d u D*»] -» 

- l -g p [ D*»t c (d p Dl d v D$ (d p Dl d„D;) t c D*» ] p£ . 



(4) 

(5) 
(6) 

(7) 

(8) 
(9) 

(10) 



With these Lagrangians we can write and evaluate the contributions of Figs. 2-4 to the total vertex function. 



III. RESULTS 

The Tr a (q) D(p) -D* (f/) vertex function T^(p 2 ) for an off-shell D is written as 

K{P 2 ) = -r a q a G{p 2 ) , 



(11) 



where G(p 2 ) is a form factor, such that the physical coupling constant is g^DD* = G{m 2 D ). We consider two kinds 
of loop corrections to this vertex, containing pion and rho intermediate states, denoted respectively by F 7T (p 2 ) and 
F p (p ). The perturbative evaluation of these functions gives rise to divergent integrals and G(p 2 ) can be determined 
only up to yet unknown renormalization constants. 

The use of standard loop integration techniques, such as dimensional regularization and MS subtraction of diver- 
gences, for all diagrams, allows one to write the form factor as 



G(p 2 ) =K + C V FAp 2 ) + C p F p (p 2 ) 



(12) 



where K, C n and C p are constants. These constants incorporate the bare couplings go and g p , the usual parameters 
associated with renormalization and, in this work, are determined by comparing the general structure of G(p 2 ) with 
the results from QCD sum rules. 

Keeping only the terms which depend on p 2 , the explicit evaluation of the diagrams given in Figs. 2-4 yields: 
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The functions I\- xyz and Ii xy entering these results are Feynman integrals, with lower labels indicating the inter- 
mediate propagating states. Upper indices represent the Lorentz tensor structure, which is realized in terms of the 
external variables q, p and p', and defined by the relationships 
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^xyz 
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xyz 
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(15) 
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(16) 



Using (13) and (14) into (12), we obtain the form factor as a function of p 2 , the D four-momentum squared. At this 
stage, it still contains three unknown parameters, which are determined by adjusting the function G(p 2 ) to the QCD 
sum rule points taken from rcf. [16]. Those results are displayed in Fig. 5, where P 2 = —p 2 , together with our best 
fit (x 2 ~ 10~ 3 ) represented by the solid line. Computing the value of G(p 2 ) at p 2 = m 2 D , we arrive at the following 
value for the coupling constant: 



gD'Drr = 17.5 ± 1.5 , 



(17) 



in very good agreement with experiment. The errors quoted come from the QCDSR points, which contain a typical 
error of ~ 10 %. In the same figure we also show the results of the fits of the QCDRS points with two mixed 
monople-dipole structures with three free parameters, namely 



G FI (p 2 ) - C 
G FII ( P 2 ) = C 1 



A 2 - P 2 
A 2 
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A! 



'D 



A\-v 2 



m 



D 
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+ C 2 



m. 



A 2 -p 2 



(18) 



(19) 



which yield xj ~ 10 -3 (dashed line) and Xii ~ 10 -2 (dash-dotted line), respectively. Inspecting Fig. 5 one learns 
that these alternative structures, reasonable as they are, diverge significantly from the loop calculation in the region 
where the D is not too off-shell, stressing the importance of a proper hadronic treatment of the form factor in that 
region. 

In this work we also consider the B*Bir vertex with an off-shell B, which can be obtained by a straightforward 
replacement of the charmed particles with those containing the quark b in the results described above. Using the 
same Lagrangians and fitting procedure with results from ref. [16], we obtain the curve shown in Fig. 6. The resulting 
coupling constant is: 



ff B . Bw =44.7±1.0 



(20) 
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in good agreement with lattice estimates. 

As far as practical applications are concerned, our numerical results for the form factors Gd* d-k{p 2 ) and Gb* Btt(p 2 ), 
in the whole range m 2 D ^ < p 2 < — 5 GeV, are very well described by the mixed monopole-dipole structure given by 
eq.(18) with the parameters 

Gd-dAp 2 ) -> C = 8.7 ; Ai = 5.1 GeV ; A 2 = 2.9 GeV , 

Gb-*Bk(p 2 ) ->C = 22.4 ; A 1 = 7.8 GeV ; A 2 = 7.1 GeV . (21) 

Our good results allow us to believe that the use of meson loops can significantly reduce the uncertainty in the 
extrapolation of form factors, computed in the space-like region by means of QCDSR, to the time-like region, with 
the corresponding increase in the reliability of predictions for coupling constants. It is worth stressing that, apart 
from the approximations described in the introduction, our procedure has no new source of errors. 

IV. CONCLUSIONS 

We have developed a new method of improving QCDSR calculations of hadronic form factors, which consists in 
matching QCDSR results, valid mainly in the deep euclidcan region, to meson loop calculations, valid when the D is 
not too off-shell. This matching is well justified from the physical point of view, since in the intermediate and large 
Q 2 regions the relevant degrees of freedom are the quarks and gluons, with non-pcrturbativc corrections taken into 
account through the QCD condensates. The opposite happens for low values of Q 2 , where sum rules calculations 
become non-reliable due to the lack of a large mass scale. At this point, the meson exchange dynamics becomes the 
most reliable tool, but it depends on unknown constants associated with the renormalization of the mesonic vertices. 
Although the exact frontier between meson dynamics and QCDSR cannot be precisely known, the success of the 
method in the instances considered here supports the view that the matching may become useful in increasing the 
predictive power of both procedures. This encourages us to reconsider our previous form factor studies. 
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FIG. 5. The D*Dn form factor. Dots: QCDSR from [16]; solid, dash and dash-dotted lines are fits obtained with eq. (12), 
(18) and (19), respectively. 
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